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The thermal decay of linear chains from a metastable state is investigated. A crossover from rigid to 
elastic decay occurs when the number of particles, the single-particle energy barrier, or the coupling 
strength between the particles is varied. In the rigid regime, the single-particle energy barrier is small 
compared to the coupling strength, and the decay occurs via a uniform saddle-point solution, with 
all degrees of freedom decaying instantly. Increasing the barrier one enters the elastic regime, where 
the decay is due to bent saddle-point configurations using the elasticity of the chain to lower their 
activation energy. Close to the rigid-to-elastic crossover, nucleation occurs at the boundaries of the 
system. However, in large systems, a second crossover from boundary to bulk nucleation can be found 
within the elastic regime, when the single-particle energy barrier is further increased. We compute the 
decay rate in the rigid and elastic regimes within the Gaussian approximation. Around the rigid-to- 
elastic crossover, the calculations are performed beyond the steepest-descent approximation. In this 
region, the prefactor exhibits a scaling property. The theoretical results are discussed in the context of 
discrete Josephson transmission lines and pancake vortex stacks that are pinned by columnar defects. 
PACS numbers: 64.60. My, 74.50.+r, 74.80.Dm, 74.60. Ge 



I. INTRODUCTION 

The decay of metastabl&jstates in systems with oncQ 
or more degrees of freedomu(DOF) has been intensively 
studied in the last decadeaJ. The crossover from rigid 
to elastic decayoQ was studied in systems with one and 
two DOF by using methods known from the auaU»is of 
the crossover from thermal to quantum decay0il3. In 
this work we consider a system with N DOF and inves- 
tigate crossovers that occur in its thermal decay from a 
metastable state while tuning an external parameter. 

A system localized in a relative minimum of a po- 
tential energy surface can escape from the trap due to 
thermal or quantum fluctuations. At high temperatures 
the decay process is purely thermal, and most probably 
occurs through the free-energy lowest-lying saddle point 
that connects two local minima. In this paper we study 
a model where the energy surface changes upon vary- 
ing an external parameter 5. Above a critical value <5*, 
the saddle point bifurcates into new lower lying ones, 
causing an enhancement of the escape rate T. In the 
steepest-descent approximation r(<5) and its derivative 
r'((5) are continuous at <5*, whereas the second deriva- 
tive r"(<5*) diverges. This behavior can be interpreted 
in terms of a second-order phase transitionta and hence is 
called crossover of second order. 

Experimental measurements concerning the decay of 
metastable states in dc superdonducting interference de- 
vices (SQUID's) were interpreted hirterms of a saddle- 
point splitting of the potential energyO. This device con- 
sists of a superconducting ring intercepted by two Joseph- 
son junctions (JJ's). The phase differences across the 
junctions play the role of generalized coordinates. The 
inductance of the circuit couples the two phases. By re- 



ducing the bias current I that flows through the sy^em, 
the decay of the phases changes from a rigid regimeEdiiS, 
with the two phases decaying together as if they were 
rigidly coupled, to an-elastic regime, with the phases de- 
caying independently!! 

An interesting question is how such a crossover oc- 
curs in more complex systems like in a discrete Joseph- 
son transmission line (DJTL), which is a onc-dimcnsional 
array of N parallelly coupled JJ's. Instead of two DOF, 
one would then have N coupled DOF. Another exaip* 
pie of such a system is a stack of N pancake vorticestfl 
in a layered superconductor in the presence of columnar 
defectsEj. A vortex pinned by a columnar defect, but 
subject to a driving current flowing perpendicular to the 
magnetic field can escape from the trap by thermal acti- 
vation. The open question is then whether a transition 
from a rigid to an elastic behavior can be found in the 
vortex or the DJTL systems, and also if more crossovers 
inside the elastic regime would arise due to the different 
decay possibilities involving the large number of DOF. In 
this paper we analyze the crossover in the decay process 
due to a saddle-point bifurcation in systems with N > 2 
DOF. It turns out that for N = 3 the saddle points of the 
potential energy can still be solved exactly. For larger N 
we determine them pcrturbatively. Furthermore, we find 
that for N 3> 1 a second crossover from boundary to bulk 
nucleation can take place in the elastic regime. 

The thermal escape rate T t h — Pexp(— Ua/ksT) is de- 
termined in the rigid and elastic regimes for an arbitrary 
number of particles, by assuming an overdamped motion 
out of a weakly metastable state. Far from the saddle- 
point bifurcation, T t h is evaluated within the Gaussian 
approximation, including the pre-exponential factor P. 
Close to the crossover from rigid decay to boundary nu- 
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FIG. 1. A current biased one-dimensional array of identi- 
cal parallelly coupled Josephson junctions, also called a dis- 
crete Josephson transmission line. The relevant degrees of 
freedom are the phase differences ip n across the junctions. 



cleation, we calculate the rate beyond steepest descent 
and find that P displays a scaling property. 

The paper is organized as follows: In Sec. II we intro- 
duce the model that can be applied both to DJTL and to 
pancake vortices in layered superconductors in the pres- 
ence of a columnar defect. In Sec. Ill we determine the 
crossover from rigid decay to elastic boundary nucleation 
and the corresponding decay diagram. We show that in 
the elastic regime a second crossover from boundary to 
bulk nucleation can occur. We evaluate the saddle-point 
solutions and their activation energies in the three decay 
regimes. In Sec. IV the thermal escape rate is calculated. 
Finally, we discuss our results and draw our conclusions 
in Sec. V. 



II. MODEL 

A. Free energy 

Let us consider a system of N degrees of freedom 
u = (uo, ■ ■ ■ , itjv-i), each of them experiencing a single- 
particle potential U{u n ), and interacting with one an- 
other via spring-like nearest-neighbor interactions, 
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where k is the spring constant. We assume that all the 
particles are initially situated near a local minimum of 
the potential U. The coordinates u n measure the dis- 
tance of each particle n from this minimum. Close to the 
local minimum u n = 0, the single-particle potential can 
be approximated by a cubic parabola, 



U(u n ) = U E 



(2) 



Here S <C 1 is a tunable parameter. The constants Ub 
and R are the characteristic energy and length scales, 
respectively. At u n — R5 the single-particle potential 
has a maximum. The energy difference between the local 
minimum and the maximum is Ub = U(RS) = UbS 3 . 




FIG. 2. A stack of "pancakes" produced by a magnetic 
field B applied perpendicular to the layers. The pancakes are 
coupled to each other via magnetic interaction and Josephson 
currents. A columnar defect pins the vortex. When a current 
j is flowing through the system, a Lorentz force fx acts on 
the pancakes, reducing the energy barrier the vortex has to 
overcome to escape from the defect. 

Among the physical systems that can be described by 
Eq. (0) are the DJTL, see Fig. |l]. The potential energy 
of a system of N identical JJ's in the presence of a bias 
current I is 
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Here the phase differences across the JJ's are given by 
(p , . . . , ipff—i- The first term in Eq. (||) represents the 
interaction energy due to the inductances between the 
loops. Here only the self-inductances of the loops are 
taken into-account, whereas the mutual inductances are 
neglectecO. The elastic constant is K = Ej/LI^, where 
Ej = ($o/27r)/ c is the Josephson energy, I c the crit- 
ical current of a single junction, L is the inductance, 
and <f>o = hc/2e is the flux quantum. The second term 
represents the tilted washboard potentials of the driven 
JJ's that arise due to the relation between currents and 
gauge invariant phases across the junctions. If we con- 
centrate on the experimentally most interesting limit 
of currents I close to criticality, NI C — I -C NI C , the 
tilted washboard potential can be well approximated by 
its cubic expansion, and we can identify £ = V with 
U n = i fn + RS/2 - tt/2, U B = <iV2Ej/3, R = 2^2, and 

S = y/(l- I/NI C ). 

Another physical realization of the model described by 
Eq. (Q) is a stack of pancake vortices trapped in a colum- 
nar defect, which is artificially introduced in a layered 
superconductor. Both the magnetic field that produces 
the pancake vortices and the columnar defect are perpen- 
dicular to the superconducting layers; see Fig. ||. Once a 
bias current j = je y flows through the layers perpendic- 
ular to the magnetic field pointing in the z direction, the 
pancakes will be driven by the resulting Lorentz force. 
The corresponding free energy reads 
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The displacement of the nth pancake vortex from its equi- 
librium position in the columnar defect is now given by a 
two-dimensional vector u„ = (u ntX , u niy ). The first sum 
in Eq. (^) models the magnetic and Josephson couplings 
between the layers by elastic interactions between pan- 
cakes in adjacent layersE-J Here Ei = (e /7 2 ) \n(X a b / £ab) 
is the elastic constant, £o = 3>o/ (47rA {,) 2 is the vortex 
self energy, 7 = X c /X a b is the anisotropy ratio of the 
penetration depths A c and X a b, s is the inter layer spac- 
ing, and £ a b is the in-plane coherence length. The second 
sum contains the columnar defect pinning potentials U p 
felt by the single pancakes and the Lorentz force density 
fi = $0 j A e z /c, where e z is the unit vector pointing 
perpendicular to the planes. The potential U p is smooth 
on the length scale £ a {, with a local minimum at the cen- 
ter of the defect. An upper estimate for the depth of the 
potential well is given by Ub ~ feo-Lp(-R/£ ;,), where R 
is the radius of the columnar defectlij. The parameter t 
denotes the superconducting layer thickness. In the large 
current limit, 5 = yl — j/jc <C 1 gives a measure of how 
close the current j is to the critical current j c . Then 
the sum of the pinning and the Lorentz part of the free 
energy is approximately 
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where we have kept only the terms that are of order <5 3 . 
The terms proportional to S(u niV /R) 2 and v? n y u nyX /R 3 , 
that are of the order 5 4 , have been neglected. Hence 
the displacements in the y direction are essentially de- 
coupled from the displacements in the x direction. As 
a consequence, two identical integrals over u nyV appear 
in the enumerator and in the denominator of the decay 
rate expression,!] which will cancel each other. For this 
reason, we will neglect u n ,v in the following. Renaming 



we obtain Eq. (hi) with k = ei/s. 



B. Decay rate 

Well above the crossover temperature To that separates 
the thermally activated decay regime from the quantum 
tunneling regime, T ^S> Tq, the escape of the DOF from 
the pinning potential can be described by a Langevin 
equation r\ 11 +V£(u) = f(t), assuming that the motion 
is overdamped. Here 77 denotes the friction coefficient. If 
we consider the resistively shunted model for the DJTL, 
•q is the inverse shunting resistance. For the vortex_prob- 
lem, 7] is given by the Bardeen-Stephen coefficientta. 

The white noise random force f(t) represents a heat 
bath at temperature T. It has ensemble averages 



(/<(«)) = and {Ji{t)fj{t')) = 2r 1 k B T8 ij 5{t - t'). In the 
limit of weak metastability, where the barrier is much 
larger than the thermal energy U a S> fcsT, the corre- 
sponding Klein-Kramers, equation can be reduced to a 
Smoluchowsky equationEI. The escape ratCpTt/j for the 
(quasi)stationary case was determined to be 
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where u' 6 S, u 6 V, S is the hypersurface in the config- 
uration space intersecting the saddle point (s) perpendic- 
ular to the unstable direction(s), V is the configuration 
volume occupied by all metastable solutions and /ig is the 
curvature of the energy surface £(u) along the unstable 
direction evaluated at the saddle point. 

Solving Eq. in the steepest-descent approximation, 
one can derive the Arrhcnius law 



T th = P(S) exp - 



Ua(S) 

k R T 
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The activation energy U a is obtained by evaluating the 
energy functional (|l|) at the sadd le-point configuration, 



which will be done in Sec. III. The computation of 
the prefactor P is a more involving task. In this case, 
we have to analyze the spectrum of the curvature ma- 
trix d n d m £ at the minimum and at the saddle point, 
since P describes the contributions to the rate that stem 
from the fluctuations around the extrema. At a char- 
acteristic value 5 = 6* the saddle bifurcates indicating 
a crossover from a rigid regime to an elastic regime. In 
the crossover region, the steepest-descent approximation 
cannot be applied. However, even beyond the steepest- 
descent approximation, the form of Eq. (Q) remains valid. 
The calculations of P will be performed in Sec. [V. 



III. 



SADDLE-POINT SOLUTIONS AND THEIR 
ACTIVATION ENERGIES 



The thermally activated escape from the local mini- 
mum u m i n — (0, . . . , 0) of the potential proceeds mainly 
via the saddle-point solutions u s of (|l|). These unsta- 
ble stationary solutions satisfy V«£ (u s ) = 0, and their 
curvature matrix H(u s ) with elements 
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has at least one negative eigenvalue. 

A. Saddle-point bifurcation 

The saddle point u rs = (R5, . . . , RS), which we call 
the rigid saddle point (rs), can be readily identified. In 
Appendix A we calculate the eigenvalues of a curvature 
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FIG. 3. The decay diagram of a system with a small 
number of degrees of freedom TV. The solid line and the dots 
indicate the crossover from rigid to elastic decay at 8 = 8, as 
a function of TV. 



matrix for a uniform extremal solution. Using Eq. ( A10| ) 
we find the eigenvalues for H(u rs ), 
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The lowest eigenvalue [Iq s 
that there is at least one unstable direction, 
only one, if 8 is smaller than 
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However, when S — > 5*, the eigenvalue y\ s = 6Ub(S* ~ 
S)/R 2 vanishes. At S = 6* the saddle splits indicating 
the existence of an elastic saddle-point configuration u es . 
Below, we will show that for 8 > 8* the energy £(u es ) is 
smaller than £(u rs ) = NUb- Hence the elastic saddle- 
point configuration u es instead of the rigid one is the 
most probable configuration that leads to decay. One 
identifies the energy of the most probable configuration 
with the activation energy U a . The saddle-point bifurca- 
tion can thus be interpreted as a crossover between two 
types of decay: the crossover from a rigid regime with 
an activation energy U a {8 < <5*) = NUb to an elastic 
regime with U a {8 > 8*) = £(u es ). The corresponding 
decay diagram is shown in Fig. [| 



B. Rigid and elastic saddles 

We now calculate the elastic saddle-point solutions. 
First, we discuss the appearance of the elastic saddle 
in the crossover regime for arbitrarily many DOF. The 
evolution of the elastic saddle point with increasing 8 is 
elucidated by analyzing the exactly solvable case of three 



DOF. Far from the crossover, the three-particle result 
is used to make an ansatz for the TV-particle solutions, 
which can again be determined perturbatively. 

Near the crossover, we expand the elastic solution 
around the rigid one, u es = u rs + Au. Then £ is most 
conveniently represented in the coordinate system of the 
principal axis of H(u rs ), where £ is diagonal in the co- 
ordinates up to second order. The transformation is 
achieved by rewriting Au as a trigonometric polynomial, 



u n = RS 
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Here the coordinates qk are the dimensionless amplitudes 
of the Fourier modes with a wave number k that mea- 
sure the deviations from the rigid saddle-point solution 
it J*/ = R8. In this coordinate system, the energy func- 
tional reads 
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where we define qj, = for k < or k > TV, and 
q = (go, ■ • • i Qn-i)- In the new coordinate system, the 
dimensionless eigenvalues flk of the curvature matrix are 
given by p, = (R 2 /U B S)no and fi k = (R 2 /2U B S)Hk for 
k ^ 0. The different prefactors are due to transformation 
(|ll|). At the rigid saddle one finds 
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where Sq^ is the Kronecker delta function. For 5 < £*, 
where the saddle-point solution is the rigid one with u n — 
RS, all the values q n = 0. The second order expansion of 
£ around the rigid saddle point reads 
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At the crossover, jl\ s vanishes and the quadratic ap- 
proximation of £ becomes independent of q\. Since 
large fluctuations in q\ would not contribute to the free 
energy, this approximation becomes insufficient within 
the crossover regime where p,™ <C 1. Thus, in or- 
der to describe the free energy contributions of fluctu- 
ations in qi more properly, higher-order terms in q\ that 
arise due to the coupling to the other fluctuation coordi- 
nates have to be taken into account. One estimates that 
A 2 £ ~ q n ^\q m ^\ ~ q n =^iq\- In comparison, the third- 
order terms qkq m qn with to, 1 are much smaller and 
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hence can be neglected. Since q\ is only coupled to q n 
and q2, one finds 



£(q) = NU B 
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In the following, wc define the small parameter e = (1 — 
8*/ 8) — —ft™/3, which measures the distance from the 
crossover. It is positive in the elastic regime and negative 
in the rigid one. Within the crossover regime, — 1 <C 
e <C 1. By solving V£ = 0, one finds the extrema. In 
addition to the extrema already found in the rigid regime, 
an elastic saddle-point solution q es with a single kink 
emerges slightly below the crossover, for 8 > 8*, 
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AD ft™ 
= 0, k > 2, 

where ft,™, ft™, and 

D = -18/(2£5') - 9/{Aft™) = 3/2 - 9/{Aft™ 



(17) 



are evaluated at the crossover. This elastic solution has 
a lower activation energy U^ s « NUb8 3 [1 — Cc 2 (8)] than 
the stiff solution. Here C = (54 - 81/ ft™)/32D 2 is a 
positive constant of the order of unity, since ft™ > 6 for 
N > 3. Since both U a (8) and its derivative U' a (8) are 
continuous, but U"(8) is discontinuous at 8 = <5», the 
crossover from rigid to elastic decay is of second order. 

In order to illustrate that in our discrete model, close to 
the crossover, boundary nucleation is the dominant pro- 
cess leading to decay in the elastic regime, we will study a 
chain consisting of three particles, where the saddle-point 
solutions can be determined exactly. The parameter e 
can now take any value in the interval — oo < e <C 1 — 8*. 
After substituting ft™ = —6, ft™ = — 3c and ft™ = 6 — 9e, 
the free-energy function reads 
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From the extremal condition V£ = we calculate the 
extrema and find that slightly below the crossover in the 
elastic regime, only q es = (gg s , qf s , g| s ), with 
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FIG. 4. The saddle-point solutions u n of a system with 
three degrees of freedom as a function of the barrier param- 
eter (5/5*. For <5 < (5* the system escapes rigidly from the 
local minimum of the potential via a configuration where all 
the particles are sitting on top of the barrier, u n = u rs . At 
(5 = 5* the saddle splits and the elastic regime is entered for 
(5 > (5* . With increasing S, uq — Uq s approaches the minimum, 
ui = ul s — 7?<5* and the last particle U2 = 112 s — > R(S + <5») is 
hanging over the maximum of the single-particle potential. 

is a possible elastic saddle-point solution. Energetically, 
the sign in front of qf s does not have any relevance since 
qi appears only quadratically in £ . It arises due to the ex- 
istence of two degenerate solutions that can mapped into 
each other by changing the sign of q\ , which is equivalent 
to a mirror symmetry transformation. Inserting the so- 
lutions for the elastic saddle q es into £, we can represent 
the free energy as a function of e; 



£(q es ) = L> B (3-3e 2 + e 3 ). 



(22) 



At e = one finds £(q es ) = £{<lrs)- Fot c > 0, the 
value of £(q es ) is smaller than that of £(q rs ). Thus there 
is a smooth crossover from the rigid q rs to the elastic 
configuration q es , which becomes the most probable one. 
To summarize, the activation energy of a three particle 
chain is given by 
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in the rigid and elastic regimes, respectively. In order to 
visualize the most probable configuration leading to de- 
cay, we represent the saddle-point solution in the original 
coordinates uq,u\ and u-i as a function of the parameter 
8. Wc find that for 8 > <5* 
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FIG. 5. (a) Rigid saddle-point solutions, (b) Elastic sad- 
dle-point solutions. 
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Note that there exists a second solution with the same 
energy, which can be found by simply exchanging the 
indices and 2. The results are displayed in Fig. [| and 
illustrated in Fig. |^. By increasing the barrier parameter 
6 above 5* , the symmetry along the defect is broken as the 
elastic saddle-point solution develops. When S is raised 
further, particle approaches the potential minimum at 
u-min — 0. Particle 1 tries to adjust between its neighbors. 
It is dragged toward the minimum by particle 0, but, 
due to the coupling to particle 2, there will be a finite 
distance between the particles 1 and 0. On the other 
hand, particle 2 has swapped to the other side of the 
maximum. 

Far in the elastic regime, ^S> N 2 , we can generalize 
this picture to arbitrary N. Making the ansatz un-i ^> 
un -2 3> ~ we find the approximate solutions of 

V£ = 0, 
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and the equivalent saddle u n 
tion energy 
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The activation energy U a is displayed in Fig. ^ for N = 
2, 3, and 4. Note that in this limit the elasticity term 
kR 2 <C 2UbS and the activation energy resembles that of 
a single particle U a ~ Ub3 3 with a renormalized barrier 
parameter. This means that for large S the system cannot 
gain much energy by nucleating at the boundary and 
bulk excitations become important. The bulk saddles 
are particle like excitations at position m with a double 
kink, 
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where \m — n\ > 1 . They have an activation energy 
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FIG. 6. The activation energy U a normalized to the ac- 
tivation energy of a single particle UbS 3 as a function of the 
barrier parameter 8 for various number of particles N. For 
iV = 2, 3 the results are exact, for TV > 3 the activation energy 
is calculated perturbatively in the crossover regime 5 ~ 8 t and 
in the limit of large 5. The activation energy for N = 2 and 
the experimental data (full dots) are taken from Ref. 5. 



which is larger than the activation energy of the elastic 
boundary saddles. Though energetically not preferable, 
for N 3> 1 the decay can occur via bulk saddle-point so- 
lutions if the barrier parameter exceeds a crossover value 
6 > 5 0S . The crossover to this new regime will be dis- 
cussed in more detail in Sec. [V. 



IV. PREFACTOR 

Having determined the activation energies U a {6) for 
the different regimes, the remaining task is to calculate 
the prefactor P(S) in Eq. (Q). Rewritten in terms of 
q = (q , q N -i), Eq. (§) reads 
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Here q = (qq, q\, . . . , qn-i) is running along S and 
Q = (^0:1 ) is probing V. In the denominator, q < 
ensures that the integration is only performed over sta- 
ble configurations. The additional prefactor arises when 
transforming the integrals to the q system and taking 
into account that = UbS^Lq/R 2 . 

A. Far from the crossover: Gaussian approximation 

In the Gaussian approximation, the integrals in the 
numerator and in the denominator in Eq. ( |35| ) are eval- 
uated by taking into account only the quadratic fluctua- 
tions around the saddle point q s , 
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and the local minimum q„ 
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respectively. Thus one obtains a prefactor 
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where the sum over the saddle index s takes into account 
the contributions of equivalent saddles. Here (/i™ n )£t™ m 
and (An)/"n are the (dimensionless) eigenvalues of the 
curvature matrices H(u m i n ) and H(u s ) evaluated at the 
local minimum (q m m)u m m and the saddles (q s )u s , re- 
spectively. In contrast to a system with translational 
invariance, in the finite systems considered here there is 
no Goldstone mode of the critical nucleus. Hence, well 
above and below the crossover, where /if ^ 0, the evalu- 
ation of P is not corrupted by divergences. 

In the rigid regime, we take only the energetically 
lowest-lying saddle into account, and the sum over s re- 
duces to a single contribution. With the det ermin ant s 
det-ff(u m , n ) and detiJ(u rs ) given in Eqs. (A7) and (AS) 
in Appendix [A], we find 



P(8 < 6*) 



W B 8 
irrjR 2 



lh(ATfi) tanh(fi/2) 



sin(iVft) tan(Q/2) 



1/2 



(39) 



where = 2 arcsinh(tj/2) and Cl = 2 arcsin((D/2) with 
lu = uj = y/ 6U b 8/kR 2 . Below the crossover, two equiv- 
alent low-energy saddle-point solutions arise, as was dis- 
cussed in Sec. III. The sum over both saddles gives rise 
to the factor 2 in 



P(S > <5») = 2 



I Mo" I 



detiT(u min ) 



\detH(u e 



1/2 



(40) 



In Eqs. (All) and (A12) we have estimated the determi- 
nant det_ff(u es ) and the eigenvalue /z§ s , respectively, in 
the limit 8 ^> <5*. We obtain 



P(6 > 5*) 



QUbS 
nr/R 2 



[1 + 0(6*/8)] . 



(41) 



As already mentioned in Sec. Ill, for N ^> 1 a crossover 
to a regime can occur, where the decay dominantly occurs 
via bulk excitations. The number of DOF Nb s , where 
the crossover from boundary to bulk nucleation occurs, 
is found by comparing the corresponding rates according 
to Eq. (|7|). In the bulk regime, one has approximately A" 




FIG. 7. Scaling property of the prefactor P as a function 
of 8 near 8*. P/P s is shown as a function of the distance e 
from the crossover. Both P s and e s are system specific scaling 
variables. 



equivalent saddles and thus with Eqs. ( A13 ) and ( A14 ) 
the prefactor is given by 



(42) 



Comparing the rates for boundary and bulk nucleation 
with U a given by Eqs. ( |30| ) and and P given by Eqs. 
( ]4l| ) and (]4^), respectively, we obtain 



2k B T\n(N/2) 



kR 2 



1/2 



(43) 



Note, that within our approximations the choice of the 
system specific parameters N, R, k and the temperature 
T is restricted to values that meet the constraint <5f, s <C 1 . 



B. Near the saddle-point bifurcation: Beyond 
steepest descent 

In the crossover regime, where 8 — > <5* and hence 
|e| — * 0, the prefactor calculated in the Gaussian ap- 
proximation diverges as P ~ l/\/e due to the vanishing 
eigenvalue fi\ = — 3e. The divergence can be regularized 
by taking into account the third order terms in q\ in the 
approximation of £ (q') around the saddle point in Eq. 
d35|). Defining the system-dependent scaling variables 

[54tanh(Q/2) sinh(iVf!)] 1/2 U 7 J 4 S 9 J 4 



7T 3/2 71 R 3 (Nk B TD) 1 / i t&n{Tr/2N) 



and 



_ / 16k B TD 
" S ~ \9NSIU B 

we show in Appendix |b] that 



1/2 
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P(e) = P s F(e/e s ), (44) 
where the function F is found to be 

tt^M exp(j, 2 ) [l_ 1/4 (y 2 ) - 7 1/4 (2/ 2 )] , <S < o*, 

8-V*r(l/4), <5 = ^, 



exp(y 2 ) IAy 2 )+F(y 2 ) 



5 > 



(45) 



For large |e/e s | the prefactor given in Eq. (|44| ) matches 
with the Gaussian result. However, in the crossover 
regime, where |e/e J< 1, the Gaussian prefactor deviates 
strongly from Eq. (E3) , as expected, since here the Gaus- 
sian approximation becomes invalid. Since we considered 
a metastable situation, where ksT -C UbS 3 , we have 
6 S « 1. Hence, the crossover regime is extremely narrow, 
| S— 5*| <C <5*. The function F = P/P s , which is shown in 
Fig. [?], reflects two interesting aspects. First, one realizes 
that the behavior of the rate is smooth at the crossover. 
The divergences that occur in the Gaussian approxima- 
tion are regularized by taking into account higher orders 
of the fluctuation coordinates. Second, F can be regarded 
as a scaling function, where the constants e s and P s con- 
tain the system-specific parameters. The scaling rela- 
tion is universal in the sense that it does not depend on 
the details of the considered system. Of course, a con- 
straint is that the crossover must be of second order to 
guarantee the validity of the perturbative treatment that 
we applied. However, we have excluded systems with a 
single-particle potential that enforce a first-order transi- 
tion from the beginning. Note, that Eq. ([fZJ) was found 
by taking into account only the cubic terms of the modes 
qo, <7i, and qi. These long- wavelength excitations deter- 
mine the decay process at the crossover, where the dis- 
creteness of the system becomes irrelevant. Hence the re- 
sult can be applied to continuous systems as well. In fact, 
a similar crossover function is found at the second-order 
transition from thermal to ouantum decay of a single par- 
ticle in a metastable stateE2l. Formally, this theory can 
also be used to describe a rigid-to-elastic crossover in the 
thermal decay of an elastic line escaping from a homoge- 
neous defect, but with 'periodic instead of open boundary 
conditions, which we considered here. Note that the scal- 
ing function found in Ref. 19 is different from ours. One 
can indeed show, that the functional form of the scaling 
function is influenced by the symmetry of the system. 



V. DISCUSSIONS AND CONCLUSIONS 

We studied the thermal decay of a chain of elastically 
coupled particles from a metastable state. The metasta- 
bility arises from each of the particles being trapped in 
a local minimum of their single-particle potential. The 



energy barrier that separates the local minimum from en- 
ergetically lower-lying ones can be tuned by a barrier pa- 
rameter 5. At 8 = the energy barrier vanishes and the 
metastability ceases to exist. With increasing 6, we find 
three regimes. For small S, the decay occurs mainly via 
a rigid configuration, where all the DOF leave the trap 
at once. At = 2kR 2 sin 2 (tt/2N)/3Ub a saddle-point 
bifurcation occurs, which marks a crossover from rigid 
to elastic motion. For S > <5* the decay occurs mainly 
via boundary nucleation. However, at even higher values 
1 5 > Sb s > 5* a crossover to bulk nucleation can take 
place. 

Our main goal was to evaluate the thermal decay rate 
T t h = P exp(— t/ a /fcsT) in the three regimes. This in- 
volves the calculation of the prefactor P and the ac- 
tivation energy U a . The latter is given by the energy 
£ of the most probable configuration leading to decay, 
namely, the lowest-lying saddle-point solution. We solved 
the problem for N = 3 particles exactly. Furthermore, 
we treated the case of an arbitrary number N of DOF 
perturbatively in the crossover regime and deep in the 
elastic regime. We have shown how the system uses 
its elasticity to lower the activation energy in the elas- 
tic regime. Whereas in the rigid regime the activation 
barrier is U™ = NUbS 3 , in the elastic regime near the 
crossover U^ s w U™(1 — Ce 2 ), where e = 1 — 6*/ 5 and 
C ~ 1 is a positive constant that depends on the de- 
tails of the potential. Increasing 6 in the clastic regime, 
the particles first escape via nucleation at the bound- 
aries with an activation energy U^ s ~ UbS 3, + kR 2 5 2 /2, 
where the first term arises from the potential energy of 
the activated particle and the second term is the elastic 
energy of the kink that occurs in the boundary saddle. 
Due to the imposed free (von Neumann) boundary con- 
ditions, this kind of activation is energetically preferred 
compared to bulk nucleation with an activation energy 
U^ s ~ UbS 3 + kR 2 5 2 . Since the bulk saddle consists of 
two kinks, twice the elastic energy is needed to activate 
a bulk nucleation process. However, in large systems, 
with N ^> 1, bulk nucleation becomes more probable for 
1 > S > S bs = ^2k B T\n{N/2)/nR 2 . Above S bs the 
many possibilities to excite a particle somewhere in the 
bulk, which grow as N in the prefactor P, outnumber the 
two possibilities of boundary nucleation. At large 5, the 
elastic interaction between the particles becomes less and 
less important and the activation energy approaches the 
energy UbS 3 which is needed to excite a single particle 
over the barrier independently of the others. To discuss 
the relevant energy scales, we now fix all variables except 
N. The crossover occurs when the number of DOF is in- 
creased above N bs = 2exp(K_R 2 <5 2 /2fcsT). Hence, when 
the elastic coupling is weak and the temperature is high, 
bulk nucleation already occurs at lower values of N bs - 
The crossover is thus determined by the ratio of elastic 
energy and thermal energy. 

Second, we determined the prefactor P. Far from the 
rigid-to-elastic crossover, the calculation of the prefactor 
P was done in Gaussian approximation both in the rigid 
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and elastic regimes. Near the crossover, the Gaussian 
approximation breaks down due to a diverging integral, 
which is caused by a vanishing eigenvalue of the curva- 
ture matrix. By taking into account higher orders in the 
fluctuation coordinates, we remove the divergence and 
obtain a smooth behavior of the rate at the crossover. 
The prefactor of the rate exhibits a scaling property 
P/P s = F(e/e s ). The function F is universal, but de- 
pends on the symmetries of the model. The scaling pa- 
rameters P s and e s are system-specific constants. 

At the saddle-point bifurcation U a (5), U' a (S), P(S), 
P'(S), and P"(S) are continuous, whereas U£ (5) is dis- 
continuous. Hence T t h(8) and T' th (d) are continuous, but 
T" h (S) is discontinuous. Interpreting U a as a thermo- 
dynamic potential, one easily sees the analogy between 
the crossover described here and a second-order phase 
transition. This analogy becomes even clearer when the 
integral in the enumerator in Eq. (0) is interpreted as 
the reduced partition sum over the DOF transverse to 
the unstable direction. Note that close to the crossover 
the discrete structure of the model becomes unimportant, 
this kind-pf-erossover can also be also found in continuous 
systemsc2l~E3. The question arises whether first-order-like 
transitions could occur also in the thermal decay of elas- 
tic chain systems- As in the crossover from thermal to 
quantum decayoE3 the type of the crossover depends cru- 
cially on the shape of the single-particle potential U(u n ). 
For a cubic parabola as is discussed in this work, the 
crossover is of second order. However, one could imagine 
other physical systems where the single-particle potential 
has a form that causes a first order transition. 

The discrete model that we have used here is quite 
general. In the following we will discuss the application 
of the theory to two physical situations, the dynamics 
of the phases in DJTL's and the thermal creep of pan- 
cake vortices in layered superconductors with columnar 
defects. 

DJTL are parallelly coupled 

one dimensional Josephson-junction arrays, and the N 
DOF in this case are the phase differences across each 
of the N Josephson junctions. In current driven DJTL, 
metastable states occur when the DOF are trapped in 
a local minimum of the tilted washboard potential com- 
mon to these systems. For N = 2, the problem reduces to 
the decay of the phases in a current biased dc SQUIDlT B. 
Both the rigid decayed, where rfche two phases behave as a 
single one, and the elastic caseS, where the two phases de- 
cay one after another, were experimentally observed. In 
the continuous limit, N — > oo, the system becomes ideate, 
cal to a long J J. The rigid-to-elastic crossover occursE3E3 
when the junction length Lj becomes of the order of the 
Josephson length Lj ~ ttXj. Here we analyzed a model 
for a DJTL, that provides a system to study the interme- 
diate case of decay from a metastable state with a finite 
number of DOF. An experimental investigation of the 
rigid-to-elastic crossover requires that the current I can 
be driven through the crossover current 7* = NI C (1 — <5 2 ). 
An orientation for the choice of the system parameters 



can be obtained by comparison with the dc- 
noting that 7* - NI C cx h 2 c 2 / (e 2 L? I 2 N 4 ) . A system- 
atic experimental study of the rigid-to-elastic crossover 
as a function of the system parameters L, I c , and N is 
still lacking and would be highly desirable. A remaining 
question was, if additional crossovers occur in systems 
with a large number of DOF. In addition to the rigid- 
to-elastic crossover due to a saddle-point bifurcation of 
the potential energy, we find that in systems with large 
N a second crossover from boundary to bulk nucleation 
can take place. DJTL's with a large number of DOF of- 
fer the possibility to observe such a crossover by varying 
system-specific parameters or the temperature. 

Let us now discuss our theory in the context of a sin- 
gle stack of pancake vortices trapped in a columnar de- 
fect in a layered superconductor. In the presence of a 
current density j that flows within the layers, the vor- 
tices are driven by the resulting Lorentz force. Once 
thermally activated from the defect, the pancake stack 
starts to move through the sample until it is trapped by 
another defect. The resulting motion is called thermal 
vortex creep. A typical example for a layered system 
is a high-temperature superconductor (HTSC). A HTSC 
like YBCO is characterized by an anisotropy 7^5, 
and the ratio of the penetration depth to the coherence 
length is \ab/iab ~ 100. The distance between the lay- 
ers and their thickness are s ~ t ~ £ a b, and the de- 
fect radius is R ~ 2£ &. In order to observe the transi- 
tion from rigid to elastic decay experimentally, the ratio 
(ic — j*)/jc > must be sufficiently large. However, 
substituting the defect energy Ub ~ teo hi(-R/£ab) and 
the elastic energy £iR 2 /s, with e ; = (e /j 2 ) ln(A Qb /£ Qb ), 
into (J c — j*)/j c — <5», one finds that even in sys- 
tems with low anisotropy and a small number of lay- 
ers (j c — j*)/j c < 10~ 2 , indicating that the phenomenon 
could hardly be observed experimentally in high-T c su- 
perconductors since j* is very close to j c . Thus, for large 
currents j c — j ' <C j c a s considered here, the vortex sys- 
tem turns out to be mainly in the elastic regime where the 
layered structure of the material is important. Then, the 
activation barrier U a is of the order of the single-particle 
barrier Ub(1 — j/jc) 3 ^ 2 , which can be interpreted as a 
vortex creep induced by the e scap e of individual pan- 
cakes from the columnar defectoO. This "decoupling" 
regime can be also entered from the low-current half- 
loop regime j <C j c , when the width of the bulk critical 
nucleus becomes of the order of the layer separationll3. 
We find that at low temperatures T the thermal creep is 
induced by boundary (surface) nucleation. It would be 
interesting to investigate experimentally if the crossover 
from bulk to surface nucleation might be observed in thin 
layered samples. In sum, we calculated analytically the 
creep rate for coupled particles trapped in a metastable 
state and found that an interesting behavior arises from 
the interplay between elasticity, pinning, discreteness and 
finite-size effects. 
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APPENDIX A: DETERMINANT AND 
EIGENVALUES OF THE CURVATURE MATRIX 

1. Recurrence relation for the Hessian matrix 



As was shown in Sec. [V, the prefactor P of the ther- 
mal decay rate is a function of the determinant and the 
eigenvalues of the curvature matrix evaluated at the rela- 
tive minimum and the saddle points, respectively, see Eq. 
(|38|). The curvature or Hessian matrix Hjv with matrix 
elements H nm (u Q ) — d n d m £(u ) determines the nature 
of £ at the extremum Uo If all eigenvalues of Hjv(uo) 
are negative (positive), Uo is a relative maximum (mini- 
mum). If some of the eigenvalues are positive and some 
are negative, then Uo is a saddle point. For £ (uo) with 
N > 3, the Hessian matrix reads 



/flg£(u) — K 



Hjv(u) 







-an 




\ 





\ —an 











In the case of open boundary conditions a — 0, the diag- 
onal elements are given by 



d 2 n £(u) 



K + U"(u n ), n = 0,N-l 
2K+U"(u n ), < n < N- 1. 



In the discussion that follows, we introduce 



/I 



Dn = det 

N 



- X -1 

-1 2 + xi 



V o 



2 + X A r_ 2 -1 

-1 1 + XN-l J 



(Al) 



which is used to calculate both the determinant and the 
characteristic polynomial of Hjy. For example, in order 
to calculate the determinant of the normalized Hessian 
Hn/k for N > 4, one sets x n — U" {u n )J n. Below, we 



will derive a recurrence relation, which is used to deter- 
mine Dfq in some special cases. 

By shifting the last column to the first and then lifting 
the bottom row to the top, one can rewrite the determi- 
nant as 



/ 1 + xjv-i 
l + x -1 
-1 



Dn = det 

N 







-1\ 



V 





-1 











.4 



N-2 



J 



where the (N — 2) x (N — 2) matrix A^^2 is given by 
f2 + xx -1 ••• \ 

-1 2 + x 2 -1 ... : 

'•• '-. '■■ 



.4 



N-2 — 



V o 



-1 2 + 2^-3 -1 

-1 2 + x N - 2 J 



In the following, we will consider the case where x\ = 
. . . = xn-2 = x. Note that xo and xn-i can be arbitrary. 
Expanding 13 at, we find with G n = det^4„ 

D N = (1 + x N -i) [(1 + x )G N -2 - G N - 3 ] 

-(l + xo)G N - 3 + G N - 4 . (A2) 

Expanding the determinant G n according to_the last 
row of A n , one finds the recursive relatione^ G n = 
(2 + x n )G n -i~G n -2 that can be rewritten as a difference 
equation 

(G n - G„-i) - (G„_i - G„_ 2 ) - x n G n -i = 0. (A3) 

The initial conditions are given by the determinants G\ 
and G2, 



Gi = 2 + x, 

G 2 = {2 + x) 2 -1. 



(A4) 



For 2 < N < 4, we can use the recurrence relations for 
G n , if we define Go = 1, G_i = 0, and G_2 = — 1. 



2. Uniform case 

The solution of these difference equations is possible 
for special cases. We now analyze the uniform case where 
x = xq = • • • = a^jv-i- Then Eq. (A2) simplifies to 



D N = (1 + x) 2 G N -2 - 2(1 + x)G N ^ 3 + G N _ 4 . (A5) 



10 



a. Determinant at the relative minimum, x > 



3. Nonuniform case 



We first discuss the case of the local minimum u = 
Umin, where x = uj 2 > 0. Imposing the initial conditions 
given by Eq. (A4), one obtains a solutionE3 of Eq. (A3), 



T N — l 



where 



Using Eqs. 



N 



sinh(iVfi) 
sinh n 



(A6) 



sinh 

(ED and (L 

lo 2 G^"_\ = 2tanh 



2 2 
we obtain 

n 



sinh(iVfi). (A7) 



Approximate solutions for the determinant and the 
eigenvalues can be obtained deep in the elastic regime, 
(5/5* > 1. 



a. Elastic boundary saddle (d~b s > 8 S> <5„) 

For the elastic boundary saddle-point configura- 
tions obtained in Eqs. (p7j)-(p9|), to highest order in 
(5/(5* one finds that J7"(u ) = ••• = U"(un~z) ~ 
6U B S/R 2 , U"(u N _ 2 ) ~ 6U b 5/R 2 -2k, and U"{u N ^x) rs 
-W B 8/R 2 - 2k. 

With x n = U"{u n )/ k one obtains for the determinant 
up to O (S N - 2 ) 



D 



N ~(l + x N ^)(2 + x N - 2 )(l + xo)G%? 3 . (All) 



b. Determinant at the rigid saddle, x < 

In the same way as for the local minimum, one obtains 
Dn at the rigid saddle u = u rs but now with negative 
x = —Co 2 < 0. One finds 



N-l 



sin(M2) 
sin Q 



where 



and hence 



D 



. fl UJ 

sin — = — . 

2 2' 



CI 



N 



2 tan - I sin(JVn) 



(A8) 



(A9) 



The ratio D 7 ^ m /D t ^, which is needed to calculate the 
prefactor in the elastic regime is found to be 



N 



D 



= -1 



N 



kR 2 



0[{8,/5) 2 



To calculate the eigenvalues, we set again x n = 
U"(u n )/n — [i Ik. The characteristic polynomial D^{p) 
is now up to 0(S N ~ 2 ), given by 

D N (n) sb (1 + x N -i)(2 + x N - 2 ){\ + x )G N - 3 (n). 

Thus, to lowest order in (5, we find that the smallest eigen- 
value is 



Mo ~ = -« - 



6U B S 
~R 2 ~' 



(A12) 



b. Elastic bulk saddle (6 > 5b s ) 



c. Eigenvalues 

The eigenvalues of Hat are found by evaluating the 
roots of the characteristic polynomial, det(Hjy — ul) = 0. 



We have again a determinant of the form of Eq. (Al), but 
now with x n = U"(u n )/K — (a/k, such th at w e can define 



D N (ix) = k~ n det(Hjv-/iI). Using Eq. (A9) we find that 



the roots where Z?Ar(/i) = are given by fl m = mir/N, 
where m — 0, . . . , N — 1. Inserting fi m into Eq. ( |A8| ) 
yields uj m — 2sin(f2 m /2), hence D^{ii m ) — for 



= 4. sin 2 (^)+U"(u ), 



(A10) 



which are the eigenvalues of Hat(uo) for a given uniform 
extremal solution uq = (uq, . . . , uq). 



For the elastic bulk saddle-point configurations ob- 
tained in Eqs. ( |3l| ) and ( |32"| ) to highest order in 5/6* 
one finds for a double kink situated at m, U"(u m ) ~ 
~6U B 8/R 2 - 4k, U"{u m±1 ) w 6U B S/R 2 - 4k, and for 
\n-m\>l U"(u n ) « 6U B 5/R 2 . With x n = U"{u n )/n 
and using periodic boundary conditions, the determinant 
is approximately given by 

D^n ~ (2 + x m +i){2 + x rn )(2 + a; m _x)Gjy!! l 3. 



The ratio D™ in /D>£ is 



N 



AkR 2 
3U B 6 



(A13) 



The characteristic polynomial Dn([i) is now up to 
0{5 N - 2 ) given by 
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D N {n) w (2 + x m+ i)(2 + x m )(2 + x m -i)GN-a(p), 

where x n = U"(u n )/n — \xjn. Thus, to lowest order in 8, 
we find that the smallest eigenvalue is 



= -2k- 



6U B S 



APPENDIX B: PREFACTOR IN THE 
CROSSOVER REGIME 

1. Rigid regime (5<8,) 

For 8 —> 8*, both the eigenvalue /z" and the determi- 
nant D r ^ vanish. Hence the Gaussian integral containing 
H\ s in Eq. (|3^) diverges, and third-order terms in qi have 
to be taken into account. In the rigid regime, the third- 
order expansion of E in qi is given by Eq. ([l5|). The con- 
tributions to P of all degrees of freedom except qi € S 
are found by Gaussian integration: 




In the following we first derive an approximate expression 
for Yl fin" 1 / Iln^i A™ and then evaluate the remaining 
integral over qi . 

For the calculation of the product term we use the rela- 
tion n^n^/Iln^iMn = W ^n™ I ■ Let us analyze 
D r j$ for close to zero. Recall that 



Cj 2 = - 



= 4 sin — — . 

K \2NJ K 



Inserting this expression into Eq. flAS| ) in the limit of 
small ijl\ s , we find 



Mi 



AT 2Ksin( 7 r/A r ) : 
such that, to lowest order in /i\ s , 



sin(A^) 



and 




2k sin(7r/iV) : 



tan (2Ar) 



Hence 



~rs 
Ml 



— cos 2 f ^r^l tanh ( ~ ) sinh(iVfi). 



A' 



TV 



\2NJ 



The integration over qi yields 



dqi exp 



NU 



2k B T 



(A14) = L/EL 




exp 



r,rs\2 



nu b {K 

16k B TD 



K 



1/4 



~rs\2 



16k B TD 



(B2) 



where £) as defined above in Eq. (|17|) arises during the 
Gaussian integrations over qo and qi- K\u is the mod- 
ified Bessel function. We make the substitution fi\ s = 
— 3e. After defining 



riV-l ~, 



1/2 



P, 



and 



[54tanh(^/2) sinh(ATfi)] 1/2 U 1 J 4 8l /i 
n 3 / 2 rjR 3 (Nk B TD) 1 / A tan(7r/2AT) 



/l6fc s TL>\ 1/2 
\ 9A^|Z7b / 



(B3) 



(B4) 



which are constants to leading order in e, we obtain the 
prefactor of the rate for the rigid region of the crossover 
regime <5<<5* :, 



J i/4 ( ^ 




7 -i/4 72 



exp [ - 



(B5) 



2. Elastic regime (S>S t ) 

In the elastic regime near the crossover, where e>0, we 
expand £ fq) around the perturbative elastic saddle-point 
solution (fl6|), 

£ (q) = S (q es ) + ^ (2) ({&» + ^ (3) ({&}), 

where £( 2 ) and f^ 3 ' contain the terms of second and third 
order, respectively, and = Qk — are the fluctua- 
tions around the elastic saddle point. By introducing the 
shifted fluctuation coordinates for m ^ 1, 



LL rs 
rm 



with ^4o = — 3, ^2 = —3/2 and Aj>2 = 0, we find, for 
the quadratic part to leading order in e, 

m^l 

Note that p,™ are the dimensionless eigenvalues evalu- 
ated at the rigid saddle-point configuration. Within the 
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crossover regime, to leading order in e, the eigenvalues 
at the elastic saddle-point solution = p^L.^ are in- 
dependent of e, except jl\ s — —2p\ s = 2e/3. The higher- 
order contributions to the expansion read 



G 

Transforming the fluctuation coordinates a second time, 
we find 

2 



ll = 6 



-£ 2 



f (q) - £(q es ) + + T 



m^l 



By using (qf s ) 2 = -/if/21) = 3e/2.D_ we evaluate the 
integrals as in the previous paragraph]^ 



d^i exp 



£i 2 - (9i s ) 2 





AT* 


{^-1/4 


r (Ar) 2 " 


+ h/4 


r (^r) 2 " 


} 


2V2 V 




16k B TD 


\6k B TD 





x exp 



16k B TD 



where 2x/4 an d I -1/4 are modified Bessel functions. The 
prefactor of the rate for the elastic regime (5><5» in the 
crossover region then reads 
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